Ground state fidelity in bond-alternative Ising chains with Dzyaloshinskii-Moriya 

interactions 
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A systematic analysis is performed for quantum phase transitions in a bond-alternative one- 
dimensional Ising model with a Dzyaloshinskii-Moriya (DM) interaction by using the fidelity of 
ground state wave functions based on the infinite matrix product states algorithm. For an anti- 
ferromagnetic phase, the fidelity per lattice site exhibits a bifurcation, which shows spontaneous 
symmetry breaking in the system. A critical DM interaction is inversely proportional to an alter- 
nating exchange coupling strength for a quantum phase transition. Further, a finite-entanglement 
scaling of von Neumann entropy with respect to truncation dimensions gives a central charge c ~ 0.5 
at the critical point. 

PACS numbers: 75.10.Pq, 64.70.Tg, 75.30.Kz 



I. INTRODUCTION 

Recent advanced material technologies have made it 
possible to access low-dimensional quantum systems. 
Furthermore, material synthesis has offered a great op- 
portunity to explore more intriguing lower-dimensional 
spin systems rather than well-understood conventional 
spin systems In such a low-dimensional system, for 
instance, alternating bond interactions and/or less sym- 
metry interactions in spin lattices can be realizable in 
synthesizing two different magnetic atoms. Of particular 
importance, therefore, is understanding quantum phase 
transitions in which one-dimensional spin systems are un- 
likely found naturally. 

Normally, quantum fluctuations in a low-dimensional 
spin system are stronger than higher dimensional spin 
systems Q. Quantum phase transitions driven by 
stronger quantum fluctuations then exhibit more inter- 
esting and novel quantum phenomena in low-dimensional 
spin systems. The effects of alternating bond interac- 
tions, especially, have been intensively studied theoreti- 
cally in spin systems such as antiferromagnetic Heisen- 
berg chains [31, Heisenberg chains with next-nearest- 
neighbor bond alternations l9l . llOll . a tetrameric Heisen- 
berg antiferromagnetic chain and two- leg spin lad- 
ders [13, A recent experiment has demonstrated a 
realization of a bond-alternating chain by applying mag- 
netic fields in a spin- 1/2 chain antiferromagnet |14| . 

In this study, we will consider one-dimensional Ising- 
type spin chains with an alternating exchange coupling. 
Actually, this bond alternation cannot destroy the anti- 
ferromagnetic phase of the uniform bond case but just 
quantitatively changes the ground state properties orig- 
inating from a dimerization of the spin lattice. Then, a 
less symmetric interaction can play a significant role to 



induce a quantum phase transition. To see a quantum 
phase transition, we will employ a Dzyaloshinskii-Moriya 
(DM) interaction [l^ which results from the spin-orbit 
coupling. 

Based on t he g round state fidelity Tg'] with the IMPS 
presentation [17|, we discuss the quantum criticality in 
the system. It is shown that a uniform DM interaction 
can destroy the antiferromagnetic phase, which is a con- 
tinuous quantum phase transition, and its critical value 
is inversely proportional to the alternating exchange cou- 
pHng strength. 



II. MODEL AND NUMERICAL METHOD 

Let us start with a spin-1/2 Ising chain with antisym- 
metric anisotropic, and alternative bond interactions on 
the infinite-size lattice. Our system can be described by 
the spin Hamiltonian 



H — JiS^S^^i + Di ■ [Si X Si+i), 



(1) 



'Electronic address: 
t Electronic address: 



libo.phyQcqu.edu.cnl 
sycho@cqu.edu.cnl 



where Si = {Sf , Sf , Sf) are the spin operators acting 
on the i-th site. The exchange interaction is chosen as 
Ji = 1 — ( — 1)V and the alternative bond interaction 
is characterized by the relative strength r of exchange 
coupling for the even and odd lattice sites. To describe 
an antisymmetric anisotropic exchange coupling between 
the two spins on the lattice, we employ a uniform DM 
interaction Di = D, that is characterized by the DM 
vector D = {D.^,Dy, D,). For r = and Z3 = 0, Eq. (P) 
is reduced to the conventional Ising chain Hamiltonian. If 
r — and D = Dz, Eq. ([T|) can be mapped onto the XXZ 
spin chain model which has a quantum phase transition 
from the gapped Neel or antiferromagnetic (AFM) phase 
to the gapless Luttinger Liquid (LL) phase at the critical 
point Dc = 1 [l^]. This study will then be focused on 
the antiferromagnetic exchange interaction Ji > 0, i.e.. 
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^ r ^ 1, and a transverse DM interaction denoting 
(0,0,7^). 

The Hamiltonian in Eq. ([T]) is actually invariant under 
the transformation ?7 = H C^2j ^ U2i+i with U2i — for 
2i-th site and U^i+i = for (2i+ l)-th site. Our model 
Hamiltonian then possesses a Z2 symmetry generated by 
the transformation \J . The ground state of the system 
may undergo a spontaneous Z2 symmetry breaking which 
gives rise to a quantum phase transition between an or- 
dered phase and a disordered phase. 

For a quantum spin system with a finite N lat- 
tice site, its wave function with the periodic bound- 
ary condition can be expressed in the matrix prod- 
uct state (MPS) representation [3] as j^E') = 
Tr . . . Al^l] IsIilsPl . . . s[^l), where is a site- 

dependent X X X matrix with the truncation dimension 
X of the local Hilbert space at the i-th site, |sW) is a 
basis of the local Hilbert space at the z-th site, and the 
physical index s takes value 1, • • • , d with the dimension 
d of the local Hilbert space. This MPS representation 
for a finite lattice system can be extended to describe an 
infinite lattice system. To do this, for an infinite lattice, 
one may replace the matrix with fWaI'I [17], where 
pl^l is a three- index tensor and A^*' is a diagonal matrix at 
the i-th site, which is called the canonical infinite matrix 
product state (IMPS) representation. 

If system Hamiltonian is translational invariant for an 
infinite lattice, for instance, our system Hamiltonian de- 
scribe by Eq. ([T]) has a two-site translational invariance, 
the two-site translational invariance allows us to reex- 
press the Hamiltonian &s H — Eift.I*''"''^!, where is 
the nearest-neighbor two-body Hamiltonian density. In 
such a case, one can introduce a two-site translational 
invariant IMPS representation, i.e., for the even (odd) 
sites A (B), only two three-index tensors rA(S) and two 
diagonal matrices \a{B) can be applied in representing a 
system wave function: 
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FIG. 1: (Color online) Fidelity per lattice site d{D, D' = 
0.26) as a function of D for various truncation dimensions 
X with r = 0.5. For higher truncation dimension, a bifur- 
cation point Db(x) moves to saturate to its critical value. 
Inset: Extrapolation of bifurcation point Db(x)- From the 
numerical fitting function Db(x) ~ Db{oo) + bx~'^ with 
b — 0.335 and c — 1.994, the critical point is estimated as 
Dc = -Dfl(oo) = 0.264. 



erators acting on two successive i and i + 1 sites through 
the Suzuki- Trotter decomposition p^ . After absorbing 
a time-slice evolution gate, in order to recover the IMPS 
representation, a singular value decomposition of a ma- 
trix is performed, which contracted from Ta, ^b, one 
Xa , two Xb and the evolution operators. Then only 
the X largest singular values are retained. This proce- 
dure yields the new tensors F^, F^ and Xa that are used 
to update the tensors for all the sites. As a result, the 
translational invariance under two-site shifts is recovered. 
Repeating the above procedure until the ground-state en- 
ergy converges yields the ground-state wave function of 
the system in the IMPS representation. 



I*) 



{.W} 



■ FaA^FbAbFaAaFbAb ■ 



, sHs['+i]s[»+2] 



.[^+3] 



(2) 

Note that, actually, for an infinite lattice sites, the diago- 
nal elements of the matrix Xi are the normalized Schmidt 
decomposition coefficients of the bipartition between the 
semi-infinite chains L(— 00, i) and R{i + 1, cxd). 

In order to find a ground state of our system in the 
IMPS representation, the infinite time-evolving block 
decimation (iTEBD) algorithm introduced by Vidal [13 
is employed. For a given initial state l'I'(O)) and 
the Hamiltonian H, a ground-state wave function can 
be yield by the imaginary time evolution |^'(t)) — 
exp[-iJr] |vI'(0))/| exp(-iJr)|*(0))| for a large enough 
T. To realize the imaginary time evolution operation nu- 
merically, the imaginary time r is divided into the time 
slices 6t = t/N and a sequence of the time slice evo- 
lution gates approximates the continuous time evolution. 
Meanwhile, the time evolution operator exp [— /I'^'^+^ljr] 
for (5r ^ 1 is expanded to a product of the evolution op- 



hi; 



FIDELITY PER LATTICE SITE AND 
BIFURCATIONS 



One can define a fidelity F{D,D') = |(*(D')I*(^))I 
from a ground state wavefunction '^[D). A fidelity per 
lattice site (FLS) d [16] can be defined as 



\nd{D,D') = lim 



\nF{D,D') 
N ■ 



(3) 



where N is the system size. Remarkably, FLS is well 
defined and plays a similar role to an order parameter 
although F{D, D') trivially becomes zero in the thermo- 
dynamic limit when N approaches infinity. The FLS sat- 
isfies the properties inherited from fidelity F{D,D'): (i) 
normalization d{D,D) = 1, (ii) symmetry d{D,D') = 
d{D',D), and (iii) range < d{D,D') < 1. 

On adapting the transfer matrix approach pl| . the 
IMPS representation of the ground-state wave-functions 
allows us to calculate the fidelity per lattice site (FLS) 
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d(D,D'). Let us choose as a reference state for 

the FLS d{D,D'). For D' = 0.26, in Fig. [H the FLS 
d(D, D') is plotted as a function of the DM interaction 
parameter D for various values of the truncation dimen- 
sions with a randomly chosen initial state in the IMPS 
representation. Figure [1] shows a singular behavior of 
the FLS ^(I?, 0.26), which indicates that there occurs a 
quantum phase transition across the singular point. A 
bifurcation behavior of the FLS 0.26) is also shown 
when the interaction parameter D becomes smaller than 
its characteristic singular value that can be called a 'bi- 
furcation point' Db- The bifurcation points depend on 
the truncation dimension x, i.e., Db — Db{x)- the 
truncation dimension x increases, the bifurcation occurs 
starting at the lower value of D. For x — >■ cx), the bi- 
furcation point Db{oo) at which a bifurcation starts to 
occur can be extrapolated. In the inset of Fig. [U we use 
an extrapolation function Db{x) — o, + bx~'^i character- 
ized by the coefficients a, h, and c being a positive real 
number, which guarantees that Db{oo) becomes a finite 
value. The numerical fitting gives a = 0.264, b = 0.335, 
and c = 1.994. From the extrapolation, the bifurcation 
point is shown to saturate to Db (oo) ~ a which can be re- 
garded as a critical point Dc = Db (oo) Actually, 
the critical point also corresponds to the pinch point of 
the FLS in the thermodynamics limit, i.e., x ~^ oo. Con- 
sequently, a FLS bifurcation point Db{x) plays the role 
of a pseudo phase transition point for a given finite trun- 
cation dimension x in the MPS representation. In addi- 
tion, the continuous function behavior of the FLS across 
the bifurcation point implies that a continuous quantum 
phase transition occurs at the critical point [l6| . 



In Fig. [U the bifurcation occurs for D < Db{x): 
which is captured in the IMPS representation with a 
randomly chosen initial state. In fact, the bifurcation 
behavior of the FLS means that there are two possi- 
ble ground states for D < Db{x) while there is a sin- 
gle ground state for D > Db{x)- Such a property of 
the ground states can be understood by the Z2 symme- 
try of the Hamiltonian from the invariant transformation 
UHU^ = H . In the thermodynamic limit, there are two 
possible ground states satisfying UHU'' = H, that is, VPg 
or U'i'g. Once a spontaneous symmetry breaking hap- 
pens, the system can choose one of two possible ground 
states or U'^g, which indicate a broken symmetry 
phase. For symmetric phase, the system has a single 
ground state, being a linear combination of two possible 
ground state, which should satisfy the transformation in- 
variance. In Fig.[Tl then, the FLS is plotted from the two 
fidelities, i.e., F = |(^'g(0.26)|^'g(D))| (upper lines) and 
F= |(*g(0.26)|J7|*g(i:'))| (lower lines). Then, the bifur- 
cation point is the transition point between the symmetry 
phase and the broken-symmetry phase. 



IV. PHASE DIAGRAM 

From the IMPS with the numerical extrapolation of the 
bifurcation points, in Fig. [5J we draw the ground-state 
phase diagram in the interaction parameter (r, D) plane. 
Below the phase boundary (red solid line) the system is 
in an antiferromagnetic phase while above the boundary 
the system is in a disordered phase. The phase diagram 
shows that the Dc is inversely proportional to the r^.. A 
best fitting function (dotted line) of the critical points 
(re, Dc) is given by Dc ~ {a{a + l))/(v^ + a) - a with 
a single parameter a = 3.4. The characteristic phase 
boundary can be understood as follows: If r ^ and 
D = 0, the lattice sites are dimerised and the Hamilto- 
nian in Eq. ([ij has a two-site translational invariance 
due to the alternating bond coupling r. As assumed, 
for the antiferromagnetic exchange interaction > 0, 
i.e., for < r < 1, the Ising chain with the alternat- 
ing bond coupling is in an antiferromagnetic state even 
though the lattice sites are strongly dimerised. If Z? 7^ 0, 
the antisymmetric anisotropic DM interaction can de- 
stroy the antiferromagnetic order originating from the 
antiferromagnetic exchange interaction Jj. The antifer- 
romagnetic order may be destroyed more easily by the 
uniform antisymmetric anisotropic DM interaction for 
the dimerised lattice sites than for the Ising chain with- 
out the alternating bond coupling because the antiferro- 
magnetic correlation between the sites becomes weaker 
due to the dimerised lattice sites. Then, to destroy the 
antiferromagnetic order, a stronger dimerisation of the 
lattice sites (bigger r) requires a much weaker uniform 
antisymmetric anisotropic interaction (much smaller D). 
Consequently, the phase boundary separating the anti- 
ferromagnetic phase and a disordered phase might have 
a inversely proportional relation between Dc and re- 

Very recently, a real space renormalization group 
(RSRG) approach [23] has been applied in the same 
model and has shown that a symmetric phase bound- 
ary is given by Dc — \/\ — (blue dashed line in Fig. 
[5]). Compared with the RSRG approach, our results from 
the IMPS show that there is a quite significant discrep- 
ancy in the phase boundary line because, contrast to the 
IMPS, the RSRG is an approximate method based on 
low-energy states, which implies that it does not cap- 
ture properly a contribution from relevant higher energy 
states. Then, the ground state phase diagram from our 
fidelity approach based on the IMPS is more reliable and 
accurate. 



V. QUANTUM ENTANGLEMENT AND PHASE 
TRANSITION 

In order to understand more clearly the quantum phase 
transition in our system, let us consider quantum entan- 
glement that can also detect a quantum phase transition 
|25| . To quantify the quantum entanglement, we em- 
ploy the von Neumann entropy, which is a good measure 
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FIG. 2: (Color online) Ground state phase diagram in 
the plane of the DM interaction and the alternating bond 
strengths. From the IMPS, the phase boundary shows that 
the Dc is inversely proportional to the ^/rZ- A best fitting 
is denoted by the dotted line. This is in sharp contrast with 
the results from the real space renormolaztion group (RSRG) 
method HI]. 

of bipartite entanglement between two subsystems of a 
pure state [26], because our ground states are in a pure 
state. Then, the spin chain can be partitioned into the 
two parts denoted by the left semi-infinite chain L and 
the right semi-infinite chain R. The von Neumann en- 
tropy is defined as S* = — Trg^ log2 Ql = ^ Trpfl log2 Qb. 
in terms of the reduced density matrix ql or pn of the 
subsystems L and R. In the IMPS representation, the 
von Neumann entropy for the semi-infinite chains L or R 
becomes 

^. = -EA2„log2A2„, (4) 

Q = l 

where A^^a's are diagonal elements of the matrix A that 
could be directly obtained in the IMPS algorithm. This is 
because, when one partitions the two semi-infinite chains 
L(— oo, •••,«) and R(i -I- 1, • • • , oo), one gets the Schmidt 
decomposition j^*) = X]a=i \,a\4'L)\4'R) ■ From the spec- 
tral decomposition, Af ^ are actually eigenvalues of the 
reduced density matrices for the two semi-infinite chains 
L and R. In our two-site translational invariant IMPS 
representation, there are two Schmidt coefficient matri- 
ces A^ and As that describe two possible ways of the 
partitions, i.e., one is on the odd sites, the other is on 
the even sites. From the A^i and As, one can obtain the 
two von Neumann entropies depending on the odd- or 
even-site partitions. 

In Fig. [3l we plot the von Neumann entropy 5' as a 
function of the DM interaction strength D for even (2i 
site) and odd {(2i + 1) site) bonds with r — 0.5. The 
von Neumann entropy for the odd bonds is always larger 
than those for the even bonds because, by the definition, 
the odd-site exchange interaction is stronger that the 
even-site exchange interaction J2i+i- Furthermore, it is 
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FIG. 3: (Color online) The von Neumann entropy S between 
left and right halves of a chain as a function of D for r = 0.5 
and X = 32. Both the von Neumann entropies for odd and 
even bonds show a singularity at Dc = 0.264. 



shown that both the entropies for the even and odd bonds 
have a singularity at the same value of the DM interaction 
strength D. Note that the singularities of the entropies 
occur at the critical point Dc — 0.264. This result then 
shows clearly that both the FLS d and the von Neumann 
entropy S give the same phase transition point. As a 
consequence, in fact, the von Neumann entropy S gives 
the same phase diagram from the FLS in Fig. [5] 

As discussed, for the antiferromagnetic state of our 
system, there are two possible ground states that are 
connected by the unitary transformation [/. From the 
bifurcation of FLS, then, one might expect a bifurcation 
in the von Neumann entropy too. However, contrary to 
the FLS d, in Fig. [Sj no bifurcation is seen in the von 
Neumann entropy S for the antiferromagnetic state even 
though the initial state is randomly chosen in the IMPS 
representation. The reason for the absence of bifurcation 
in the von Neumann entropy is why the singular values 
Ai in Eq. ^ do not depend on the unitary transforma- 
tion because the unitary transformation U acts only on a 
single site of the spin lattice in the IMPS representation. 



VI. CENTRAL CHARGE AND UNIVERSALITY 
CLASS 

At a critical point, characteristic singular behaviors of 
thermodynamics system properties depend only on few 
features such as dimensionality and symmetry, which can 
be classified by the concept of universality classes. Espe- 
cially, the central charge can be used for the classification 
of universality classes 27, 28]. Owing to implement the 
IMPS representation, we can obtain a central charge c 
and a so-called finite-entanglement scaling exponent k 
numerically via the unique behaviors of the correlation 
length ^ and the von Neumann entropy S with res pec t 
to the truncation dimension x at a critical point (29l.l30j| , 
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FIG. 4: (Color online) (a) Correlation length ^ as a function of 
the truncation dimension x s-t the critical point. The power 
curve fitting ^ = ax'' yields a = 0.112 and k = 2.132. (b) 
Scaling of the von Neumann entropy 5* with the truncation 
dimension x the critical point. For k = 2.132 from (a), the 
linear fitting S = {ck/6) log2 x + ^ yields the central charge 
c ~ 0.494. Here, the alternating bond strength is chosen as 
r = 0.5. 

i.e., 

e = ax", (5) 
S = ylog2X- (6) 

In Fig. m the correlation length ^ and the von Neu- 
mann entropy S* as a function of the truncation dimension 
X at the critical point Dc for r = 0.5. Here, the trun- 
cation dimensions are taken as x = 4, 6, 8, 12, 16, 24, 32, 
and 48. It is shown that both the correlation length ^ and 
the von Neumann entropy S diverges as the truncation 
dimension x increases. From a power-law fitting on the 



correlation length ^, we have n = 2.132 and a = 0.112. 
As shown in Fig. [?] (b), our numerical result demon- 
strates a linear scaling behavior, which gives a central 
charge c ~ 0.494 with k = 2.132. Our central charge 
is close to the exact value c — 1/2. Consequently, the 
quantum phase transition in our system is in the same 
universality class as the quantum transverse field Ising 
model. 



VII. SUMMARY 

Quantum phase transitions have been investigated in 
the Ising chain with the Dzyaloshinskii-Moriya inter- 
action as well as the alternating bond-coupling. The 
FLS and its bifurcation have clearly shown a charac- 
teristic singular point as a signature of the quantum 
phase transition and behaves as a continuous function, 
which shows a continuous phase transition occurring at 
the critical point. The phase diagram was obtained from 
the FLS and the von Neumann entropy. With a finite- 
entanglement scaling of the von Neumann entropy with 
respect to the truncation dimension in the IMPS repre- 
sentation, a central charge was estimated to be c ~ 0.5, 
which shows that the system is in the same universality 
class with the quantum transverse field Ising model. 
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